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We study off-resonant collective light scattering from ultracold atoms trapped in an optical lattice. 
Scattering from different atomic quantum states creates different quantum states of the scattered 
light, which can be distinguished by measurements of the spatial intensity distribution, quadra- 
ture variances, photon statistics, or spectral measurements. In particular, angle-resolved intensity 
measurements reflect global statistics of atoms (total number of radiating atoms) as well as local 
statistical quantities (single-site statistics even without an optical access to a single site) and pair 
correlations between different sites. As a striking example we consider scattering from transversally 
illuminated atoms into an optical cavity mode. For the Mott insulator state, similar to classical 
diffraction, the number of photons scattered into a cavity is zero due to destructive interference, 
while for the superfluid state it is nonzero and proportional to the number of atoms. Moreover, we 
demonstrate that light scattering into a standing-wave cavity has a nontrivial angle dependence, 
including the appearance of narrow features at angles, where classical diffraction predicts zero. The 
measurement procedure corresponds to the quantum non-demolition (QND) measurement of various 
atomic variables by observing light. 

PACS numbers: 03.75.Lm, 42.50.-p, 05.30.Jp, 32.80.Pj 



I. INTRODUCTION 

Ever since the first generation of Bose-Einstein conden- 
sates (BEC), it has been a central task to study quan- 
tum properties of such degenerate gases. Surprisingly, it 
turned out that many properties are well explained by the 
Gross-Pitaevskii equation, which is a an effective nonlin- 
ear single-particle equation and allows to calculate the 
evolution of the average atomic density and phase. The 
density can be observed by simple absorption images af- 
ter expansion, and the phase can be mapped onto density 
modulations in interferometric setups. The limited valid- 
ity of such mean-field descriptions became apparent with 
the advent of optical lattices [l|, Q , where one has quan- 
tum phase transitions between states of similar average 
density but radically different quantum fluctuations. 

The majority of methods to characterize quantum 
properties of degenerate gases are based on matter-wave 
interference between atoms released from a trap in time- 
of-flight measurements 0, H, 0] j which destroys the sys- 
tem. Recently, a method of "Bragg spectroscopy" based 
on stimulated scattering of matter waves by laser pulses 
was applied to homogeneous BECs [sl, IgJ and atoms in 
lattices [7, 8, 9, 10]. In this case, the measured quantities 
(e.g. structure factor), which carry information about 
density fluctuations, are also accessible via matter- wave 
interference. Although the scattered light and stimulated 
matter waves can be entangled and mutually carry the 
statistical information [ill, > laser flelds are simply 
considered as a tool to stimulate matter waves. 
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In contrast to those works, the nondestructive meth- 
ods based on measurements of light flelds only, without 
destro ying at oms, were proposed in Refs. p^. Il4. [Tsl. [l6l. 
[i3,[li,[lll23,[2l|J^ for homogeneous BECs in traps and 
optical lattices [23]. Here the average amplitude of the 
scattered light is solely determined by the average atomic 
density, while the photon number and other higher or- 
der fleld expectation values contain quantum statistical 
properties of atoms. 

In this paper, we show that this is of even greater sig- 
nificance for atoms in lattices, where different quantum 
phases show qualitatively distinct light scattering. Here 
we extend the preliminary results presented in our previ- 
ous letter, Ref. [l^]. In particular, linear scattering can 
create entangled states of light and manybody atomic 
states, exhibiting a nontrivial connection of the field am- 
plitude and intensity. As a practical consequence, we 
demonstrate the possibility to distinguish between dif- 
ferent quantum phases, e.g., Mott insulator (MI) and 
superfiuid (SF), by measuring properties of a scattered 
off-resonant beam. This possibility is exhibited in sev- 
eral different ways involving simple intensity measure- 
ments, or more involved measurements of quadrature 
variances, photon statistics, as well as phase-sensitive or 
spectral measurements. A careful analysis of the scat- 
tered light provides information about global statistics 
(related to atom number at a lattice region illuminated 
by the probe), local quantities (refiecting statistics at a 
single site even without an optical single-site access) , and 
pair correlations between different sites. 

Note that we consider off-resonant and almost nonde- 
structive light scattering, corresponding to the quantum 
non-demolition (QND) measurements of various atom- 
number functions. In principle, it can be repeatedly or 
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even continuously applied to the same sample. This is 
very different from noise spectroscopy in absorption im- 
ages [1| where observations of quantum fluctuations of 
the atomic density were recently reported. 

For homogeneous BECs [H [li, [IE [H, [13, d, E, 
[20l . [2ll . |2^ . the scattered light was shown to consist of 
two contributions: the strong classical part insensitive to 
atomic fluctuations, and weaker one, which carries infor- 
mation about atom statistics. For a large atom number, 
the classical part completely dominates the second one, 
which, in some papers, even led to a conclusion about the 
impossibility of distinguishing between different atomic 
states by intensity measurements, and, hence, to a ne- 
cessity to measure photon statistics. 

In our work, we show that light scattering from atoms 
in optical lattices has essentially different and advanta- 
geous characteristics in contrast to scattering from homo- 
geneous BECs. For example, the problem of suppressing 
the strong classical part of scattering has a natural so- 
lution: in the directions of classical diffraction minima, 
the expectation value of the light amplitude is zero, while 
the intensity (photon number) is nonzero and therefore 
directly reflects density fluctuations. Furthermore, in an 
optical lattice, the signal is sensitive not only to the pe- 
riodic density distribution, but also to the periodic dis- 
tribution of density fluctuations, giving an access to even 
very small nonlocal pair correlations, which is possible by 
measuring light in the directions of diffraction maxima. 

As free space light scattering from a small sample can 
be weak, it might be selectively enhanced by a cavity. 
The corresponding light scattering from an optical lat- 
tice exhibits a complicated angle dependence and narrow 
angle-resolved features appear at angles, where classical 
diffraction cannot exist. In experiments, such a nontrivial 
angle dependence can help in the separation between the 
signal reflecting atom statistics from a technical noise. 

Joining the paradigms of two broad fields of quantum 
physics, cavity quantum electrodynamics (QED) and ul- 
tracold gases, will enable new investigations of both light 
and matter at ultimate quantum levels, which only re- 
cently became experimentally possible [H, [2^, [23| • Here 
we predict effects accessible in such novel setups. 

Experimentally, diffraction (Bragg scattering) of light 
from classical atoms in optical lattices was considered, for 
example, in Refs. [H, 0, [s^l • In our work, we are essen- 
tially focused on the properties of scattering from ultra- 
cold lattice atoms with quantized center-of-mass motion. 

The paper is organized as follows. In Sec. II, a gen- 
eral theoretical model of light scattering from atoms in 
an optical lattice is developed taking into account atom 
tunneling between neighboring sites. In Sec. Ill, we sig- 
nificantly reduce the model to the case of a deep lattice 
and give a classical analogy of light diffraction on a quan- 
tum lattice. Section IV presents a relation between atom 
statistics and different characteristics of scattered light: 
intensity and amplitude, quadratures, photon statistics, 
and phase-sensitive and spectral characteristics. Proper- 
ties of different atomic states are summarized in Sec. V. 




FIG. 1: (Color online) Setup. Atoms in a lattice are illumi- 
nated by a pump wave at angle ^o; scattered (probe) light is 
collected by a cavity at angle 9i and measured by a detector. 

In Sec. VI, we present a simple example of the model 
developed: light scattering from a lattice in an optical 
cavity pumped orthogonally to the axis. The main results 
are discussed in Sec. VII and summarized in Sec. VIII. 



II. GENERAL MODEL 

We consider an ensemble of A'' two-level atoms in an 
optical lattice with M sites. Except the presence of a 
trapping lattice potential, atoms are illuminated by and 
scatter field modes at different directions. A possible ex- 
perimental realization is shown in Fig. [TJ Here, a lattice 
is illuminated by a "pump" beam, whereas measurements 
are carried out in one of the scattered modes, which is 
treated as a "probe." Note that different experimental 
setups are possible: the modes can be either in free space, 
or selected by traveling- or standing-wave cavities, or 
even correspond to different modes of the same cavity. 
For definiteness, we will consider the case, where mode 
functions are determined by cavities, whose axes direc- 
tions can be varied with respect to the lattice axis (the 
simplest case of two standing-wave cavities at angles 6q 
and 9i is shown in Fig. [T]). Instead of varying the an- 
gles, the mode wavelengths can be varied with respect 
to the wavelength of a trapping beam. We also assume, 
that not all M lattice sites are necessarily illuminated by 
additional modes, but some region with K < M sites. 

The manybody Hamiltonian in the second quantized 
form is given by 

H = Hf + Ha, (la) 
Hf = ^ hujiajai - {viai - iliaj), (lb) 

Ha = J £r¥{v)Hai^{r) 

^27ra^ r ^3j.^t(i.)^t(j.)^(r)*(r). (Ic) 
m J 

In the field part of the Hamiltonian Hf, ai are the an- 
nihilation operators of light modes with the frequencies 
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oji , wave vectors k; , and mode functions ui (r) , which can 
be pumped by coherent fields with amphtudes rji. In 
the atom part, Ha, ^'(r) is the atomic matter-field op- 
erator, as is the s-wave scattering length characterizing 
the direct interatomic interaction, and Hai is the atomic 
part of the single-particle Hamiltonian Hi, which in the 
rotating-wave and dipole approximation has a form 



this expansion in Eq. ([T]) with Hai (HI), we get 



M 



U 



Hal 



2ma 



huj, 



Hi^Hf+Hai,{2a) 
- ihgoYla+aiUiir) - H. c.](2b) 



where the coefficients J^j correspond to the quantum mo- 
tion of atoms in the classical potential and are typical for 
the Bose-Hubbard Hamiltonian JjI : 



Here, p and r are the momentum and position operators 
of an atom of mass and resonance frequency LUa , cr^ , 
(T~, and az are the raising, lowering, and population dif- 
ference operators, go is the atom-light coupling constant. 

We will consider essentially nonresonant interaction 
where the light-atom detunings A;a ^ lui — LUa are much 
larger than the spontaneous emission rate and Rabi fre- 
quencies goo-i- Thus, in the Heisenberg equations ob- 
tained from the single- atom Hamiltonian Hi az can 
be set to —1 (approximation of linear dipoles). Moreover, 
the polarization can be adiabatically eliminated and 
expressed via the fields a;. An effective single-particle 
Hamiltonian that gives the corresponding Heisenberg 
equation for a; can be written as ifieg = Hf + Hai with 



J^!,. = J drwir ~ rO + Ki(r)) ^v{r -r,). (5) 

However, in contrast to the usual Bose-Hubbard model, 
one has new terms depending on the coefficients j'™, 
which describe an additional contribution arising from 
the presence of light modes: 



J(7 = / drw(r - r,)u,*(r)u™(r)u;(r - rj). (6) 



^•^i - ^ + ^'^ii^) + ^9o2^ X ■ (3) 

Here, we have also added a classical trapping poten- 
tial of the lattice, Vci(r), corresponds to a strong clas- 
sical standing wave. This potential can be, of course, 
derived from one of the modes a; — ad [in this case 
Vci{r) — ?i5Q|aciUci(r)P/AciQ], and it can scatter hght 
into other modes. Nevertheless, at this point we will 
consider Vci(r) as an independent potential, which does 
not affect light scattering of other modes that can be sig- 
nificantly detuned from Oci [i.e. the interference terms 
between Oci and other modes are not considered in the 
last term of Eq. ([3])]. The later inclusion of the light scat- 
tered by the trapping wave will not constitute a difficulty, 
due to the linearity of dipoles assumed in this model. 

We will consider scattering of weak modes from the 
atoms in a deep lattice. So, the fields ai are assumed 
much weaker than the field forming the potential Vci{r). 
To derive the generalized Bose-Hubbard Hamiltonian we 
expand the field operator ^'(r) in Eq. using localized 
Wannier functions corresponding to Vci(r) and keeping 
only the lowest vibrational state at each site: ^'(r) — 
biw{r — Ti), where bi is the annihilation operator of 
an atom at the site i with the coordinate . Substituting 



In the last term of Eq. only the on-site interaction 
was taken into account and U = AirasTT? /ma J dr\w(r)\^. 

As a usual approximation, we consider atom tunneling 
being possible only to the nearest neighbor sites. Thus, 
coefficients ([S]) do not depend on the site indices {Jf\ = 
Jq and Jfij-i — J'^'), while coefficients ^ are still index- 
dependent. The Hamiltonian ^ then reads 

H = Hf + JS'iV + r'B + hgl E ^ f E 

l,m \2— 1 / 

\ ( K \ jj M 

+Woj:t^ e 47^'^. +|E"^("^-i)'(7) 

t,m \<*ij> / *— 1 

where < i,j > denotes the sum over neighboring pairs, 
fii = b\hi is the atom number operator at the i-th 
site, and B — X]f=i ^l^i+i + H.c. While the total atom 
number determined by iV = X]f=i '^j conserved, the 
atom number at the illuminated sites, determined by 
Nfc = X^^Li is '^ot necessarily a conserved quantity. 

The Heisenberg equations for a; and hi can be obtained 
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from the Hamiltonian ^ as 



K 



K 



ai 



<h3> 



K 



l.m 



-J (j^' + %o'E|^4"-i)^-i(-8b) 



In Eq. ([8a[) for the electromagnetic fields a;, two last 
terms in the parentheses correspond to the phase shift of 
the light mode due to nonresonant dispersion (the second 
term) and due to tunneling to neighboring sites (the third 
one). The second term in Eq. ((8a|) describes scattering of 
all modes into a; , while the forth term takes into account 
corrections to such scattering associated with tunneling 
due to the presence of additional light fields. In Eq. (|8b[) 
for the matter field operators bi, the first term gives the 
phase of the matter-field at the site i, the second and 
third terms describe the coupling to neighboring sites. 

It is important to underline that except for the direct 
coupling between neighboring sites, which is usual for 
the standard Bose-Hubbard model, Eqs. ^ also take 
into account long-range correlations between sites, which 
do not decrease with the distance and are provided by 
the common light modes a; that are determined by the 
whole set of matter- field operators bi. Such nonlocal cor- 
relations between the operators bi, which are introduced 
by the general Eqs. ([S]), can give rise to new many-body 
effects beyond predictions of the standard Bose-Hubbard 
model [13 . 



III. SCATTERING FROM A DEEP LATTICE 
AND CLASSICAL ANALOGY 

We will significantly reduce the general model de- 
scribed by the Hamiltonian ([7]) and Heisenberg equa- 
tions ([5]). In contrast to our paper [3]| and works on 
so-called "Bragg spectroscopy" 0, H, [9l. Il0j| , we will not 
consider excitations of the lattice by light and stimula- 
tion of matter waves. The focus of the present paper is 
a study of properties of light scattered from the atoms 
in a prescribed quantum state, which is not necessarily 



the ground one. The main result is the demonstration 
of the possibility to distinguish between different atomic 
quantum states of different statistics by measuring light 
only. 

We consider a deep lattice formed by a strong classi- 
cal potential Vci(r), so that the overlap between Wannier 
functions in Eqs. ^ and ^ is small. Thus, we can ne- 
glect the contribution of tunneling to the scattered light 
by putting J^' = and Jj™ = for i 7^ j. Under this 
approximation, the matter-wave dynamics is not essen- 
tial for light scattering. In experiments, such situation 
can be realized because the time scale of light measure- 
ments can be much faster than the time scale of atomic 
tunneling. One of the well-known advantages of the op- 
tical lattices is their extremely high tunability. Thus, 
tuning the lattice potential, tunneling can be made very 
slow [ll . On the other hand, the rate of the photon escape 
from the cavity depends on the cavity relaxation rate and 
photon number, while the letter is determined by vari- 
ous parameters as atom-field detuning, pump amplitude, 
and atom number. In modern experimental setups, all of 
these parameters, especially the ultracold atom number, 
can be tuned in a very broad range. Moreover, even with 
no tunneling, various atomic quantum states with essen- 
tially different statistics can be realized until the atoms 
will decay to the ground state. Since in this paper we 
do not require the atoms to be in a ground state, light 
scattering can refiect different atomic statistics even for 
negligible tunneling. 

In a deep lattice, the on-site coefficients j'™ ^ can be 
approximated as J^'™ = it;*(ri)u,„(ri) neglecting details 
of the atomic localization. Such details are accessible 
even from the classical consideration 28, 29, 30]. In this 
paper, we will focus on essentially quantum aspects of 
the problem. 

For simplicity, we will consider scattering of a single 
mode flo ("pump"), considered as a given operator, into 
another mode ai with the relaxation rate k included phe- 
nomenologically. The Heisenberg equation (|8ap for the 
scattered light ai then reads 



2 K 

Ala 



i=l 



ai 



* Aoa 



K 



^u\{T,)U(i{Yi)T 



KGi + 771, 



(9) 



where we do not add the Langevin noise term, since we 
will be interested in normal ordered quantities only. In 
the Heisenberg equation for the matter-field operators bi 
(I8bp . only the first term is nonzero. This term affects only 
the phase of the matter field, but not the atom number 
operators hi = b\bi. Hence, though the the matter-field 
phase still depends on the common light mode, the op- 
erators fii, appearing in Eq. ([5]), are constant in time. 

We assume that the dispersion shift of the cavity mode 
gp/Aia X]i!=i is much smaller than k or detun- 
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ing between the pump and scattered light Aqi = ujq—uji. 
Thus, a stationary solution of Eq. ^ has a form 



ai = CD, C = 



iglao 



Aoa(K - iAoi) 



K 



(10a) 
(10b) 



where we have also assumed no additional pumping {rji = 
0) and replaced the operators ao,i(i) by their slowly vary- 
ing envelopes ao,i(i) [ao.i{t) = ao,i exp(— iwo^)] skipping 
in the following notations all tilde signs. 

Expressing the light operators in terms of the atomic 
ones in Eq. pH)) is a central result here, which we will use 
to study the properties of the scattered field. The depen- 
dence of the light Heisenberg operators on the atomic op- 
erators reflects the entanglement between light and mat- 
ter during the light-matter interaction. We will assume 
the pump mode to be in the coherent state, which enables 
us to consider the quantity ao as a c-number. 

In the following, we will consider a ID lattice of the pe- 
riod d with atoms trapped at Xm — rnd (m — 1, 2, .., M). 
The result for the field operator oi (|10ap with the oper- 
ator D (|10b[) has an analogy in classical diffraction. For 
scattering of a traveling wave oq in the direction of a 
traveling wave oi from a lattice with (hi) = n at each 
site, the expectation value of the field is given by 



K 



C{D) ^CJ2 



AmSkx d I 



sin(a_/2) 



(11) 



where a. 



5krrd, and 5k„ 



(ko 



kism 6q 



sin^i) is the projection of the difference between two 
wave vectors on the lattice direction, 0o,i are the angles 
between wave vectors and a vector normal to the lattice 
direction (cf. Fig. [T]), k — to / c for ujq — uji — to. 

Equation pT|) simply describes classical diffraction of 
the traveling wave ao on a diffraction grating formed by 
equally spaced atoms with positions of diffraction max- 
ima and minima (i.e. scattering angles 9i) determined 
by the parameter a_ depending on the geometry of inci- 
dent and scattered waves and diffraction grating through 
^0, |ko.i|i and d. A more general form of the operator D 
given by Eq. (jlOb[) describes also diffraction of a standing 
wave ao into another mode ai, which can be formed, for 
example, by a standing-wave or ring optical cavity. 

Equation (fTT|) shows that the expectation value of the 
scattered field is sensitive only to the mean number of 
atoms per site n and reflects a direct analogy of light 
scattering from a classical diffraction grating. Neverthe- 
less, the photon number (intensity) and photon statistics 
of the field ai are sensitive to higher moments of the 
number operators fii as well as to the quantum corre- 
lations between different lattice sites, which determines 
quantum statistical properties of ultracold atoms in an 
optical lattice and will be considered in the next sections. 



IV. RELATION BETWEEN QUANTUM 
STATISTICS OF ATOMS AND 
CHARACTERISTICS OF SCATTERED LIGHT 

A. Probing quantum statistics by intensity and 
amplitude measurements 



According to Eq. (|10ap . the expectation value of the 
photon number a|ai is proportional to the expectation 
value of the operator D*D. We introduce coefficients 
^1(^*0, ^1) responsible for the geometry of the problem: 



K 



b^^Aifii, A,{9o,0i) = ul{ri)uo{ri), 

i=l 

K 

A{6o,e^) = Y,A,{eo,e^), (12) 

where uo,i(rm) = exp(imA;o,ia;'^ + 0o,im) for traveling 
waves, and wq. 1(1*771 

) = cos(TOfco,ia;d + 00,1m) for stand- 
ing waves (m = 1,2, ...M), fco^ia; = |ko,i| sin 00,1, are 
the angles between mode wave vectors and a vector nor- 
mal to the lattice axis; in the plane-wave approximation, 
additional phases 0o,im are m-independent. 

The expectation values of D and D*D then read 

K 

{b)=Y,A,{n,) ^nA, (13a) 

7 = 1 

K 

{b*b)= J2 4*^.(«.'b-> (13b) 

K 

= (n„nfc)|A|2 + ((n2)_(^^^,))^|A,|2, (I3c) 

7=1 

R{eo,ei) = (b*b)^\(b)\^ ^ 

K 

({hah,) - 7^2)|A|2 + ((n^) - {hah,)) ^ \A,\' (13d) 

7=1 

K 

= {6ha8h,)\Af + ((<57i2) - {ihM,)) \A^\^■ (13e) 

7=1 

In Eqs. (|13p we have used the following assumptions 
about the atomic quantum state \^)'. (i) the expecta- 
tion values of the atom number at all sites are the same, 
(n-i) = n (thus, the expectation value of atom number at 
K sites is {Nk) — Nk = nK), (ii) the nonlocal pair cor- 
relations between atom numbers at different sites {hihj) 
are equal to each other for any i ^ j and will be denoted 
as {hah,) (with a ^ b). The latter assumption is valid 
for a deep lattice. We also introduced the fluctuation 
operators Shi = fM — n, which gives {Sh'^) equal to the 
variance (An^)^ 



Equation (|13ap reflects the fact that the expectation 
value of the field amplitude (jlOap is sensitive only to the 
mean atom numbers and displays the angle dependence 
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of classical diffraction given by the factor A{9q, 6i), which 
depends on the mode angles and displays pronounced 
diffraction maxima and minima. Equation l|13bp shows 
that the number of scattered photons (intensity) at some 
angle is determined by the density-density correlations. 
In the simplest case of two traveling waves, the prefac- 
tors A*Aj = eyi^[i5kx{xj — Xi)] with 5kx = fcox — kix- In 
this case, Eq. (jl3b|) gives the so-called structure factor 
(function), which was considered in the works on light 
scattering from homogeneous BEC (TtI . [Tsj . Here we es- 
sentially focus on optical lattices. Moreover, it will be 
shown, that the more general Eq. (|13bp . which includes 
scattering of standing waves, contains new measurable 
features different from those of a usual structure factor. 

Equation p3cp shows, that the angle dependence of 
the scattered intensity consists of two contributions. The 
first term has an angle dependence 1^(6*0, 611) p identical 
to that of the expectation value of the field amplitude 
squared (|13ap . The second term is proportional to the 
quantity (n^) — {fiahb) giving quantum fluctuations and 
has a completely different angle dependence X)i=i l^jP- 
The expression (I13c[ ) has a form similar to the one con- 
sidered in papers [isl. [isl . [iGj on light scattering from 
a homogeneous BEC, where the scattered intensity con- 
sisted of two parts: "coherent" (i.e. depending on the 
average density) and "incoherent" one (i.e. depending 
on the density fluctuations). Nevertheless, in the present 
case of a periodic lattice, this similarity would be exact 
only in a particular case where there are no nonlocal pair 
correlations {fiafib) = naUb = v? {{ShaSfib) = 0), which 
in general is not true and leads to observable difference 
between states with and without pair correlations. 

Further insight into a physical role of nonlocal pair 
correlations can be obtained from Eqs. (fT3dl) and dTSe]) 
for the "noise quantity" R{9o,Oi) = {I)*D) - \{D)\'^, 
where we have subtracted the classical (averaged) con- 
tribution |(-D)P to the intensity {D*D). Equation 
shows that, in the noise quantity, a term with the classical 
angular distribution \A{9o,di)\'^ appears only if the pair 
correlations are nonzero. The physical meaning of this 
result is that, in an optical lattice, it is not only the den- 
sity distribution that displays spatial periodic structure 
leading to diffraction scattering, but also the distribution 
of number fluctuations themselves. In the framework of 
our assumption about equal pair correlations, the spatial 
distribution of fluctuations {SfiaShb) can be either the 
same as the density distribution (with a lattice period d) 
or zero. In the former case, pair correlations contribute 
to the first term in Eqs. (|I3d[) and (|13ep with classical 
distribution \A{0o,9i)\'^ , in the latter case, (SfiaShb) = 0, 
and the only signal in the noise quantity is due to on- 
site fluctuations {6n'^) with a different angle dependence 
J2i=i l^iP- Note that, in general, the spatial distribution 
of fluctuations can be different from that of the average 
density and can have a period proportional to the lattice 
period d. This will lead to additional peaks in the angu- 
lar distribution of the noise quantity (jl3dp . (|I3ep . The 



generalization of those formulas is straightforward. 

Even with spatially incoherent pump ao, the intensity 
of the scattered mode a\ai is sensitive to the on-site atom 
statistics. To model this situation, the quantum expecta- 
tion value {D*D) (I13bp should be additionally averaged 
over random phases 0o,im appearing in the definition of 
mode functions in Eq. (H^). In Eq. (jl3bp . only terms 
with i = j will then survive and the final result reads 

{D*D)i^,=PoK{n^), (14) 

where pq is equal to 1 for two traveling waves, 1/2 for 
a configuration with one standing wave, and 1/4, when 
both modes ao,i are standing waves. 

B. Quadrature measurements 

The photon number a\ai is determined by the expecta- 
tion value {D*D), whereas (D) gives the field (oi) (llOap . 
While photon numbers can be directly measured, a field 
(oi) measurement requires a homodyne scheme. Such a 
measurement then makes (D) experimentally accessible. 
Actually for a quantum field only the expectation values 
of quadratures of ai that are Hermitian operators and 
can be measured. Using Eq. (|10ap and the commutation 
relation [ai, a\] = 1, the quadrature operator and its 
variance (AX^)^ can be written as 

X^^^ (aie-** + ate'*) = \C\X^_^^, (15a) 
X| = i + |Cp(l^%J^(15b) 
{AX^f EE {XD - {X^r = i + \C\'{AX^_^J^il5c) 

where C = |C| exp(i0c') and the quadratures of D are 

X^^^[De-^^ + D*e^f'), (16a) 

{Ax^r ^ (ix^r) - {x^r. (leb) 

In Eqs. (fT5|) . the phase is related to the homodyne 
reference phase, while 0c is determined by the phase of 
the pump ao and parameters of the field-matter system 
[cf. Eq. (|10ap ]. Hence, the phase f3 = (j) — (f)c entering 
Eqs. (|15p can be controlled by varying the phase differ- 
ence between the pump and homodyne fields. 

Using Eq. (fT2|) . the quadrature operator X^ reads 

K 

X^-^Afn,, Af(0o,^i)= l^zl cos (0^.-/3), 

i=l 

K 

A^{0o,9,) = Y,A^{0o,S,),{n) 

where Ai = \Ai \ exp(z0y!i.), and we defined new quantities 
Af (00,^1) and (00,^1). 
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Since Eq. HT]) for Xj^ and Eq. for 1) have a similar 
structure, the Eqs. ^3]) for the quantities {D), {b*b), 
and i? can be rewritten for the quantities {X^') , {{Xp Y), 
and (AX^)^, respectively, with the change of parame- 
ters A, (6*0, 6*1) and ^(6*0, 6*1) to Af (6*0, 6I1) and A^(6lo, 6I1). 
Thus, the above discussion of Eqs. can be repeated 
in terms of the quadrature operators with the only dif- 
ference that coefficients A^{9q,0i) and A^{9q^9i) now 
depend also on the homodyne phase. An advantage of 
this reformulation is that the expectation value of the 
non-Hermitian operator ai , which determines {D) , is now 
replaced by the expectation value of the Hermitian opera- 
tor X^, which is consistent with a procedure of measuring 
quadratures of the quantum field ai. The well-known re- 
lations between intracavity and outcoupled fields can be 
found, e.g., in Ref. [33 | for linear systems, which is the 
case in our work. 

I 



where we assumed again that correlations do not depend 
on site indices, and sites with the indices a, 5, c, and 
d are different. In Eq. (fTO|) . each prefactor containing 
geometrical coefficients Ai (|12p determines different angle 
dependences of a corresponding term. 

Thus, varying the geometry of a problem (e.g. angles 
of two modes, wavelengths of the modes or that of trap- 
ping potential determining the lattice period), one has 
access to different statistical quantities characterizing the 
quantum state of ultracold atoms. 

D. Phase-sensitive and spectral measurements 

In the derivation of Eq. (|10p . we have neglected the 

term 5o/^iaSi=i |wi(ri)P"'j i^i Eq. ^ related to the 
refractive index of atoms for the scattered light. This 
term is normally very small at large detunings. However, 
if the scattered mode is confined in a very good opti- 
cal resonator, the light experiences a very long effective 
path within the atoms, and this term shifts the phase of 



C. Photon number fluctuations 

While the intensity of the scattered light is sensitive to 
the second moments of the number operators n^, quan- 
tum statistics of the field reflexes the forth-order mo- 
ments. The variance (Ariph)^ of the photon number 
JT-ph — a\ai is given by 

(Anph)^ = (n^h) - ("ph)^ (An^i^) : -H(nph) 

= |C|4((73*2£)2) _ {^b*bf) + \C\^{b*b), (18) 

where : (An^J (af a?) - {a\ai)'' = \C\^{{b*'' b'') - 

{b* b)"^) is a normal ordered photon-number variance. 
Thus, the problem is reduced to measurements of the 
photon number \C\'^{b*b) and quantity \C\'^ {b*"^ b"^) , 
which after straightforward calculations is given by 



the scattered light. In a steady state approximation it 
amounts to the dispersion shift of a cavity mode. 

Equation ^ shows that even in the absence of the 
pump field (oq = 0), quantum fluctuations of the atom 
number enter the phase via operators of K illuminated 
sites, which depend on the atomic quantum state. In the 
simplest case of a traveling wave, |ui(ri)| — 1, and the 

operator X)t=i ~ is the number of atoms at K 
sites. As will be discussed below, in the Mott insulator 
state, the expectation value of this quantity Nk = {Nk) 
does not fluctuate. In the superfluid state with K = M, 
Nm is equal to the total number of atoms iV, and also 
is fixed. However ior K < M, Nk fluctuates strongly 
and, as will be discussed, for K ^ M corresponds to a 
coherent state with (Nj.) = (Nk)^ + (Nk)- Those sta- 
tistical properties of the atomic states are reflected in the 
phase of the light field. In particular, measurements of 
the dispersion shift of a cavity mode will show a frequency 
distribution reflecting the distribution of atom numbers. 

This also opens an alternative spectral method of de- 
termining the quantum state of the atoms in a cavity 



K 



{b*^b'^) ^ 



K 



K 



V 2 



{nanbUcnd) + 2 



K 



K 



-I- c.c. 



{2{nanbncnd) - 3{nlnbnc) + (nlnb)) 



{-{nanbUcnd) + {n 



lnbnc)) + 2^JA,\'^ {{n^ 



nbUcUd) - 2{nlnbnc) + {nlnl)) 



K 



E^? 



{{nanbUcUd) - 2{nlnbnc) + {nlnl)) + 4 



E^^ 



K 



K 



\Ai\'^{-{nanbncnd) + {nlnbUc)) 

1=1 

\A,\'{~Q{nanbncnd) + I2{nlnbnc) ~ ^{nlnb) - i{nlnl) + (n^)), (19) 

i=i 

I 
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TABLE I: Statistical quantities of typical atomic states. 



with two degenerate modes. Let us consider the mode og 
as a dynamical quantity obeying an equation as Eq. ([9]), 
which can be obtained from the set of Eq. ([8a|l . while the 
second degenerate mode is called ai. The atoms lead to 
the collective normal-mode splitting of two cavity modes 
as recently experimentally observed [s^. If the coupling 
coefficient between two degenerate modes, which is equal 
to ^p/Aoa-D [cf. Eq. ([9]) and the definition of the oper- 
ator D in Eq. p^ ]. exceeds the cavity relaxation rate 
K, a spectral doublet instead of single maximum can be 
observed in the spectrum of the output light. 

It is quite expected that the collective strong coupling 
between the modes and thus the spectral splitting de- 
pends on the number of atoms in a lattice. Interestingly, 
from the equations for ao and ai, it can be shown, that 
parameters of the normal-mode splitting (e.g. splitting 
frequency, linewidths) also depend on the atomic quan- 
tum state. So, spectral mode-splitting measurements 
also can be used to distinguish between atomic quan- 
tum phases and allow a nondestructive measurement of a 
quantum phase transition dynamics. In the following we 
will, however, restrict our study to single frequency mea- 
surements and leave a more detailed analysis of phase- 
and frequency-sensitive phenomena to other works (s^ . 



V. QUANTUM STATISTICAL PROPERTIES OF 
TYPICAL ATOMIC DISTRIBUTIONS 

Let us briefly summarize some key statistical proper- 
ties of typical states of N atoms at M lattice sites, i.e: 
the Mott insulator state (MI), superfluid state (SF), and 
a multisitc coherent-state approximation to the SF state 
(cf. Table m. 

The MI state represents a simple product of local Fock 
states at each site with precisely Ui atoms at a site i. As 
a consequence, atom numbers at each site hi (as well as 
the number of atoms at K sites Nk) do not fluctuate, 
and there is no quantum correlations between sites. 

Similarly to the pair correlations, all two-, three-, and 
four-site quantities in Eq. p9|) factorize. From the hght- 
scattering point of view, this is the most classical atomic 
state, which corresponds to periodically ordered pointlike 
atoms. We will further consider the commensurate flUing 



with rii = N/M atoms at each site, neglecting possible 
random vacancies. This can be made if one has some ad- 
ditional information that quantum fluctuations dominate 
over other, thermal or technical, sources of noise. 

The SF state corresponds to a BEC where each atom 
is in the zero quasi-momentum Bloch-state of the lowest 
band and is equally delocalized over all sites. Hence, the 
atom numbers at a given site (and the number of atoms 
ai K < M sites) fluctuate. As a consequence of the total 
atom number conservation, the numbers of particles at 
two different sites a ^ h are anticorrelated. All two-, 
three-, and four-site quantities in Eq. also do not 
factorize. 

The expectation values in the SF state can be calcu- 
lated using normal ordering and the following relations: 



h\^SF{N,M)) = 



|*sf(A^-1,M)), 



(vPsFl^l-brivI/sF) = 



N{N -1)...{N -m + l) 



where the first equation relates SFs with N and — 1 
atoms. 

We will introduce another, coherent, quantum state, 
which is often considered as an approximation to the SF 
state, and represents a product of local coherent states 
at each site. In this approximate state, the numbers of 
particles at a given site and at any K < M sites fluctuate. 
Moreover, the total number of particles at M sites is also 
a fluctuating quantity, which is a disadvantage of this 
approximation. Similarly to the MI state, correlations 
between several different sites are absent. In the coherent 
state, one has 



h\^CoY.{N,M)) = ^l — \M.cobiN,M)), 



(M/coh|6l"6ri* 



Coh) 



N" 



Comparing properties of the SF and coherent states 
in Table [H we can state that under the approximation 
N, M oo, but flnite N/M, the coherent state is a good 
approximation for local one-site quantities and correla- 
tions between different sites. Moreover, ii K <^ M, the 
SF expectation values related to the nonlocal Nk opera- 
tor are also well approximated by corresponding quanti- 
ties in the coherent state. Nevertheless, if the number of 
sites K is of the order of M, the coherent-state approxi- 
mation fails for those quantities. 

One can prove even a more general statement for the 
functions (D*D) ^ and (D*^D^) dH]), which deter- 
mine the intensity and statistics of light and are the most 
important quantities in this work. If the number of sites 
illuminated by light, K, is much smaller than the total 
number of lattice sites M, the coherent-state is a good 
approximation for calculating characteristics of scattered 
light in the limit TV, M ^ oo, but flnite N/M. If, in op- 
posite, the number of sites interacting with light is of 
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the order of the total number of sites in the lattice, this 
approximation, in general, gives wrong results. As will 
be shown, it fails for light scattering in the directions of 
diffraction maxima. The proof of the statement is based 
on the consideration of the orders of sums in Eqs. (fT3)) . 
(fT9|) . which contain geometrical coefficients Ai and are 
proportional to the powers of K, whereas factors contain- 
ing atom fluctuations have powers of M in denominators. 

Thus, light scattering from the region of a SF optical 
lattice with K <^ M sites is equivalent to the light scat- 
tering from the atoms in the coherent state (in absolute 
values both K and M can be very large). Moreover, in 
the directions outside diffraction maxima, the coherent- 
state approximation works well even in the case where 
any number of sites is illuminated. 

In the following, discussing all states, we will use the 
notations n — N/M for the atomic "density" (expec- 
tation value of the particle number at each site) and 
Nk — KN /M — TiK for the expectation value of the par- 
ticle number at K sites. These two parameters fully char- 
acterize light scattering in the MI and coherent states, 
while all three parameters N, M , and K are necessary 
to characterize scattering in the SF phase. For definitive- 
ness, we will discuss a case with large values of TV, M, 
and K where difference between odd and even number 
of lattice sites vanishes. Nevertheless, note that physical 
problems including BECs with large atom number loaded 
into lattices with small site numbers are also of great im- 
portance [H, [3^. Results for this case, can be obtained 
from expressions of this section and Eqs. and (fTO|) . 



VI. EXAMPLE: ID OPTICAL LATTICE IN A 
TRANSVERSALLY PUMPED CAVITY 

Before considering a general angular distribution of 
scattered light, we would like to present the most striking 
prediction of our model describing the difference between 
atomic quantum states, observable by light scattering. 
Let us consider a configuration of Fig. [T] where the pump 
(traveling or standing wave) is orthogonal to the lattice 
(^0 = 0), and the scattered light is collected along the 
lattice axis {Oi = 7r/2) by a standing- or traveling- wave 
cavity. This geometry coincides with the one considered 
in the context of cavity cooling [s^, [H, 111] and lattices 
in optical cavities [3l| . Atoms are assumed to be trapped 
at each lattice site [d = A/2) at field antinodes. 

In this case, the operator D is reduced to 

SfLi which, independently on an atomic 
state, gives zero for the expectation value of the field 
amplitude proportional to {D) (here we assume even 
K). This corresponds to the classical destructive inter- 
ference between atoms separated by A/2. In contrast, 
the photon number in a cavity a\ai is proportional to 
{D*b) ^ {{n^) - {nanb))K [cf. Eq. (HSc])], which is de- 
termined by statistics of a particular state, and is equal 
to zero for the MI state and to Nk for the SF state. 



Thus, atoms in a MI state scatter no photons into a 
cavity, while a SF scatters number of photons propor- 
tional to the atom number: 

(ai)Mi = (ai)sF = 0, but 
{a\ai)Mi = 0, (4ai)sF = \C\^Nk. 

Hence, already the mean photon number provides in- 
formation about a quantum state of ultracold atoms. 

The photon number fluctuations (Ariph)^ (fT8|) are also 
different for various states. In the MI state, the variance 
(A|7^|2)2 = (^*2^2^ - {D*D)^ is zero, (A|i?|2)2^j = 0, 
whereas in the SF state, Eq. (|19p gives a very strong noise 
(A|L>|2)|p = 2N]^ (in highest order of Nk)- 

Nonlinear light-matter dynamics in a cavity can lead 
to a new self-organized phase [H, H^] where all atoms 
occupy only each second site leading to doubling of the 
lattice period, d — \. The operator D (fT2)) is then re- 
duced to X]fc=i "fc = ^K- Thus, if the final self-organized 
state is a MI with cZ = A, the photon number in a cavity 
is (a|ai)scif-org = \C\'^N]^, which is proportional to the 
atom number squared and has a superradiant character. 
This result coincides with the theory of self-organization 
with classical center-of-mass motion [401 . 



VII. RESULTS AND DISCUSSION 

In the following we will compare light scattering from 
atoms in the following states: MI, SF with all sites illu- 
minated {K — M using the notation SFa/ ), and partially 
illuminated SF under the approximation A^, M —> 00, fi- 
nite n = N/M, K <C M, which will be denoted as the 
"coherent" taking into account the equivalence proved in 
Sec. V. The results for the SF/<- state with any K can be 
obtained from the general Eqs. ((13]) and ^9]) . We will 
restrict ourselves to the case of plane waves. Distinguish- 
ing between atomic states using light modes with more 
complicated spatial profiles can be analyzed by general 
expressions of Sec. IV. 

A. Two traveling waves and discussion of essential 
physics 

For two traveling waves, which can be free-space modes 
or fixed by ring cavities, the geometrical coefficients 
are Am = exp(ima_) (a_ — kQ^dsinOo — kixdsinOi), 
and Eq. (jl3ep for the noise quantity is reduced to 

R = {Snjhbf-^^^-^^^^ + iiSn^) - {ShJn,))Kl22) 
sm (a_/2) 

where the first term has the angle dependence of classical 
diffraction pT|) . and the angle dependence in the second 
term in Eq. (jl3ep is reduced to a constant (isotropic) 
one, K. In the MI and coherent states, where pair corre- 
lations {SfiaSni,) are absent, the first term is zero. In the 
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FIG. 2: (Color online) Intensity angular distributions for two 
traveling waves, the pump is transverse to the lattice (So = 0). 
(a) Intensity of classical diffraction of coherent light (curve A), 
isotropic intensity of incoherent light scattering, Eq. (|14[) . for 
coherent atomic state (line B) and MI state (line C); (b) noise 
quantity, Eq. (|23|l . for coherent atomic state (constant value 
1, line A), SF with all sites illuminated K = M (curve B), 
and MI (constant value 0, line C); (c) the same as in (b) but 
for partially illuminated SF with K = M/2. N ^ M = 30. 



MI state, on-site density fluctuations {Sn^) are also zero 
giving the zero value of the noise quantity (j22p . while in 
the coherent state, it is the on-site fluctuations {6tt?) — n 
that give isotropic contribution to R. Thus, we have 



-Rcoh 

,2 , 



-Rmi — 0, 
nK = Nk, 



'm2 sin2(a_/2) ' 



(23a) 
(23b) 

(23c) 



It is important to note, that in the SF state (|23cp . even 
in a large optical lattice with N,M oo , very small 
pair correlations (ShaSfib) = —N/M'^ can give a signifi- 
cant angle-dependent contribution to the noise quantity, 
which occurs near a diffraction maximum (a„ = 27r/, / — 
0, 1, ..), where the geometrical factor is equal to K^, and 
if the number of the illuminated sites K is of the or- 
der of M. This demonstrates the importance of nonlocal 
correlations and invalidity of the coherent-state approx- 
imation under those conditions. Outside the diffraction 
maximum, where the geometrical factor is small, pair 
correlations do not play any role and the coherent-state 
approximation works well even for all sites illuminated. 

Figure [2] shows several angle dependences of the scat- 
tered light in the case of two traveling waves. As an 
example, in all figures, we will consider atoms at each 
lattice sites providing d — Ao,i/2. In Fig. [H^a), the an- 
gular distribution of classical diffraction | {D) p (curve A) 
is shown. In the case oi d = Xq.i/2 and the pump being 
orthogonal to the lattice (^o = 0), only the zero-order 



diffraction maxima at 6i = 0, tt are possible in the clas- 
sical picture. Corresponding noise quantities R for the 
coherent (constant lines A) and SF^^ (curves B) states are 
displayed in Figs.dfb) and[2)Jc) (in MI, the noise is zero, 
which is displayed by lines C). According to Eq. (|23p . 
the intensity fluctuations are isotropic for the coherent 
atomic state, while there is suppression of intensity noise 
under scattering from the SF. The suppression occurs in 
the regions of diffraction maxima. For all sites illumi- 
nated, K = M [cf. Fig. [2Ib)], the suppression is total, 
while for K = M/2 it is only partial [cf. Fig. ^c)]. 
Outside the maxima, the dependence for SF^f is well ap- 
proximated by that for the coherent state for any K. 

It is important to underline, that in a broad range 
of angles, the number of scattered photons from the SF 
(or coherent) state is nonzero, even if the expectation 
value of the electromagnetic field vanishes, which man- 
ifests the appearance of nonclassical entanglement be- 
tween the light and manybody atomic system. Moreover, 
in contrast to MI state, atoms in SF state scatter photons 
at angles, where the classical diffraction does not exist. 

For example, in a simple configuration considered in 
Sec. VI where the pump is orthogonal to the lattices 
{00 = 0), and the scattered light is collected by a cav- 
ity along the lattice axis (6*1 — tt/2), the atoms in the MI 
state scatter no photons as in classical diffraction min- 
imum. In contrast, atoms in the SFk state scatter the 
number of photons a\ai = \C\^{D*D) = \C\^Nk, pro- 
portional to the number of the atoms illuminated [cf. 
Eq. ^ and Fig.^b) at the angle 9i =n/2]. 

For two traveling waves, the expression for D (fT^ . in 
a diffraction maximum where all atoms radiate in phase 
with each other and a_ = 27rZ, is reduced to the op- 
erator Nk- Thus, the quantity (D) = Nk = nK is 
the expectation value of the atom number at K sites 
and proportional to the average atom number at a single 
site. The intensity of the light scattered into a diffrac- 
tion maximum is determined by {D*D) = {N]^), while 
noise R — (ANk)'^ gives the atom number variance at 
K sites. The latter statement corresponds to Figs. (Hb) 
and m^c) displaying the total noise suppression in SF m 
state, where the total atom number at all sites K = M 
does not fluctuate, while for K < M, Nk is a fluctuating 
quantity and the noise suppression is only partial. 

At the angle of a classical diffraction "minimum" (for 
K ^ 1 this is approximately valid for any angle out- 
side narrow regions of maxima), the expectation value of 
the fleld amplitude is zero, as well as the flrst terms in 
Eqs. (HSc]), (fT3d)l . (fBe)) . and both the intensity (D*b) 
and noise R are proportional to the quantity {n^) — {nafib) 
giving the difference between local and nonlocal fluctua- 
tions. For two traveling waves, the coefficient of propor- 
tionality is isotropic and equal to K [cf. Eq. (|22p ] . 

For scattering of incoherent light (|14p . the intensity is 
proportional to the local quantity (n^) and is shown in 
Fig.mja) for MI (curve C) and coherent, almost the same 
as in SF, (curve B) states. This quantity can be also ob- 
tained under coherent scattering of two traveling waves. 
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if one tunes tfie angles such that the geometrical factor of 
the first term in Eq. is equal to K. Practically, this 
variant is easy to achieve only for a diffraction pattern 
with diffraction maxima, which are not too narrow. 

Hence, in an optical experiment, both global statisti- 
cal quantities related to K < M sites, local quantities 
reflecting statistics at a single site, and pair correlations 
can be obtained. It is important, that local statistics can 
be determined by global measurements, i.e., an optical 
access to a single site is not necessary. 

Therefore, light scattering gives a possibility to dis- 
tinguish different quantum states of ultracold atoms. As 
demonstrated by Eq. (|^ and Fig. O MI and SFm states 
are distinguishable in diffraction "minima" and in inco- 
herent light, while they are indistinguishable (for trav- 
eling waves) in maxima, because the total atom number 
contributing to the maximum does not fluctuate. The 
SFjvf and coherent states can be distinguished in diffrac- 
tion maxima only. The MI and coherent states can be 
distinguished in any angle of the scattering pattern. 

Measurements of the noise quantity discussed or, alter- 
natively, related quantities for quadratures (|16p or pho- 
ton number variance ([18]) , give the values, which are dif- 
ferent in orders of the emitter number Nj^ for different 
quantum states. Nevertheless, for large Nk, there could 
be practical problems in the subtraction of large values 
in a diffraction maximum to get the noise contribution. 
In some papers, a similar problem even led to a con- 
clusion about state indistinguishability by intensity mea- 
surements in BEC [l3 . [T9I |20| and, hence, to a necessity 
to measure photon statistics. A rather involved method 
to suppress the strong classical part of scattering using a 
dark-state resonance in BEC was proposed in Ref. 16]. 
In contrast to homogeneous ensembles, in optical lattices, 
this problem has a natural solution: measurements out- 
side diffraction maxima are free of the strong classical-like 
part and thus directly reflect density fluctuations. 



B. Physical interpretation and role of the 
entanglement between light and matter 

The classical analogy of the difference in light scatter- 
ing from different atomic states consists in various den- 
sity fluctuations in different states. In particular, classi- 
cal density fluctuations would also lead to impossibility 
of obtaining a perfect diffraction minimum, where contri- 
butions from all sites should precisely cancel each other. 

Scattering at diffraction maxima can be treated as su- 
perradiant one, since the intensity of the scattered light is 
proportional to the number of phase-synchronized emit- 
ters squared N^. In diffraction minima, destructive in- 
terference leads to the total (subradiant) suppression of 
coherent radiation for MI state; whereas for SF^ state, 
the intensity is nonzero and proportional to the number 
of emitters Nk, which is analogous to the emission of 
independent (non-phase-synchronized) atoms. 

Nevertheless, the quantum treatment gives a deeper 



insight into the problem. 

The expression for the SF state in Table U can be 
rewritten in the following from: 

I*sf) = 77=777X1 \/ —r-r — rl9i'92,-'7j\/), 
(VM)^ ^ y gi!g2!---gA/! 

where the sum is taken over all qi such that g,; = N. 
It shows that the SF state is a quantum superposition of 
all possible multisite Fock states corresponding to all pos- 
sible distributions of N atoms at M lattice sites. Under 
the light-matter interaction, the Fock states correspond- 
ing to different atomic distributions become entangled to 
scattered light of different phases and amplitudes [4il.[4^. 

For example, in a simple case of two atoms at two sites, 
|*sf) = l/2|2,0) + l/%/2|l,l) + l/2|0,2). In the example 
configuration of Sec. VI, where the orthogonal pump il- 
luminates lattice sites separated by A/2 (diffraction min- 
imum), the wave function of the whole light- matter sys- 
tem reads 

|*mattcr-light) = 1/2|2, 0) I7) -|- 1/V2\1, 1) |0) 

+l/2|0,2)|-7). 

Here, if we assume that the distribution |2, 0) is entan- 
gled to the coherent state of light I7), the distribution 
|0, 2) will be entangled to the similar light state with 
the opposite phase | — 7), and the distribution |1, 1) will 
be entangled to the vacuum field |0), because the fields 
emitted by two atoms cancel each other. 

In contrast to the classical case, light fields entan- 
gled to various atomic distributions do not interfere with 
each other, which is due to the orthogonality of the Fock 
states, providing a sort of which-path information. This 
leads to a difference from the classical (or MI with the 
only Fock state |1, 1)) case and nonzero expectation value 
of the photon number even in the diffraction minimum 
(I7P/2 in the above example). The absence of inter- 
ference gives also an insight into the similarity of scat- 
tering from the SF state to the scattering from indepen- 
dent (non-phase-synchronized) atoms, where interference 
is also absent. 

We would like to mention that light measurements con- 
sidered here correspond to the QND measurements [i^ of 
atomic variables observing light. Considering the Hamil- 
tonian ([7]) with neglected tunneling, it is straightforward 
to show that for the "probe observable" (i.e. measured 
light a\ai) and "signal observable" D (jlObp . all four con- 
ditions of the QND scheme summarized in Ref. [i^ are 
fulfilled. Thus, observing light one can determine var- 
ious atomic functions corresponding to the geometry- 
dependent operator D (|10b[) in a QND way. For example, 
the total atom number or atom number at some lattice 
region can be nondestructively measured in a diffraction 
maximum, while the difference between atom numbers at 
odd and even sites can be nondestructively determined 
in a diffraction minimum. 
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C. Standing waves 

If at least one of the modes is a standing wave, the 
angle dependence of the noise becomes richer. In an ex- 
periment, this configuration corresponds to a case where 
the scattered light is collected by a standing-wave cav- 
ity, whose axis can by tuned with respect to the lat- 
tice axis [3l|. Except for the appearance of new clas- 
sical diffraction maxima represented by the first terms in 
Eqs. (|13c[) . (jl3dp . (|13ep . which depend on the phase pa- 
rameters a± = kQ^dsinOQ ± fcij-dsinf?!, the angle depen- 
dence of the second term is also not an isotropic one, as it 
was for two traveling waves. This second, "noise," term 
includes a sum of the geometrical coefficients squared, 
which is equivalent to the effective doubling of the lat- 
tice period (or doubling of the light frequency) and leads 
to the appearance of new spatial harmonics in the light 
angular distribution. Such period doubling leads to the 
appearance of the peaks in the noise distribution at the 
angles, where classical diffraction does not exists. 

In Fig. [21^ a) , angular distributions of the scattered light 
are shown for a traveling- wave pump, which is almost 
orthogonal to a lattice (6*0 = O.Itt), while the probe 
is a standing wave. Classical diffraction pattern [cf. 
Fig. [3fal)] is determined by \A\'^ through the parameters 
a± and shows zero-order diffraction maxima in transmis- 
sion {9i = 9o and its counterpart due to the presence of 
the standing-wave cavity at 9i = n + 6q) and reflection 
(^1 = IT — 9q and the counterpart at — ^o)- The inten- 
sity noise for atoms in the coherent state [cf. Fig. [2Ia2)] 
is determined by X^iLi l^jP through another parameter 
2ai = 2fci2,c?sin0i and has different characteristic fea- 
tures at 9i = 0, TT, and ±7r/2. It is the latter feature 
that corresponds to the effective frequency doubling and 
appears at an angle, where classical diffraction has a min- 
imum. In the case of SFm state [cf. Fig. [3]Ja3)], pair cor- 
relations in Eqs. (|13d[) and (|13ep are nonzero, hence, both 
geometrical factors contribute to the noise distribution, 
which has the features at angles characteristic to both 
classical scattering and the light noise of the coherent- 
state case. Outside the characteristic features, the noise 
distribution is isotropic and takes a nonzero value similar 
to the case of two traveling waves [cf. Fig. [2]. Figure[3l^b) 
shows a simpler situation, where the pump is precisely 
orthogonal to the lattices (6'o — 0). 

In Fig. [3l^c), a situation similar to Fig. [3la) is shown 
for the case where both the pump and probe are stand- 
ing waves. While classical diffraction still depends on 
the parameters a±, the factor X^iLi l^iP determining 
the intensity noise depends on four parameters 2q!o,i = 
2fco,i2;(i sin 00.1 and 2a±. Thus, in the light noise from 
a lattice in the coherent and SF states, the features are 
placed at the positions of classical zero-order diffraction 
maxima and the angles, which would correspond to the 
classical scattering from a lattice with a doubled period 
d — \, where the appearance of first-order diffraction 
maxima is possible. Similar to Fig. [31[a), features at 
6*1 = 0,7r, and 7r/2 also exist. In the case 6*0 = 0, the 
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FIG. 3: (Color online) Intensity angular distributions for scat- 
tering into a standing-wave cavity, (a) Traveling-wave pump 
at Qq = O.Itt; (b) traveling or standing-wave pump at So = 0; 
(c) standing- wave pump at So = O.Itt. Intensities of classical 
diffraction are shown in Figs, (al), (bl), and (cl); noise quan- 
tities for coherent state are shown in Figs. (a2), (b2), and (c2) 
and for SF in Figs. (a3), (bS), and (c3). iV = M = JC = 30. 
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angular distribution for two standing waves is identical 
to that of one standing wave shown in Fig. E^b). 



In the SFm state, there are two types of diffraction 
maxima. In the first one, the noise can be completely 
suppressed due to the total atom number conservation, 
similarly to the case of traveling waves. This occurs, 
if the condition of the maximum is fulfilled for both of 
two traveling waves forming a single standing wave [cf. 
Fig. [2Ib)]. In the second type, even for K = M, only 
partial noise suppression is possible, since only one of the 
traveling waves is in a maximum, while another one, be- 
ing in a minimum, produces the noise [cf. Figs.[3I^a) and 
[31[c)]. In contrast to two traveling modes, in the second 
type of maxima, one can distinguish between SFm and 
MI states, since MI produces no noise in any direction. 



-;r/2 



D. Quadratures and photon statistics 



An analysis of the angular distribution of the quadra- 
ture variance (AX^)^ (jl6bp shows, that even for two 
traveling waves, new peaks due to effective period dou- 
bling appear [see Fig.SI^a)]. Additionally, the amplitude 
of noise features can be varied by the phase difference be- 
tween the pump and homodyne beams which is shown 
in Figs.lHb) for the coherent and in Fig.UJ^c) SFm states. 
In the coherent state, all peaks are very sensitive to f3. In 
the SFm state, the noise suppression at diffraction max- 
ima is insensitive to variations of /3, whereas other peaks 
are /3-dependent. The relation of {AX^')^ to the quadra- 
ture variance of the light field is given by Eq. pScp . 



FIG. 4: (Color online) Quadrature angular distributions for 
two traveling waves, (a) Quadrature for classical difTraction 
(al), quadrature variance for coherent (a2) and SF (a3) states, 
pump-homodyne phase difference /3 = 0; (b) quadrature vari- 
ance for coherent state for /3 = n/4 (bl), /? = n/2 (b2), 
/3 = 37r/4 (b3); (c) quadrature variance for SF for f3 — n/i 
(cl), /? = n/2 (c2), /? = 37r/4 (c3). Oo = 0, N = M = K = 30. 



The angle dependence of the variance (A|Z)|^)^ — 
{D*'^!)^) — {D*D)^, which is proportional to the nor- 
mal ordered photon-number variance and determines the 
light statistics (jlSp , also shows anisotropic features due to 
frequency doubling even for two traveling waves (Fig. [5]). 
In this case, Eq. (|19p is reduced to 
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FIG. 5: (Color online) Angular distributions of photon- 
number variances for two traveling waves, (a) Intensity of 
classical diffraction (curve A), isotropic intensity of incoher- 
ent light scattering, Eq. (|14|) . for coherent atomic state (line 
B) and MI state (line C); (b) normal ordered photon- number 
variance for coherent atomic state under scattering of coher- 
ent, Eq. (|24|l (curve A), and incoherent (line B) light; (c) 
normal ordered photon-number variance for SF state under 
scattering of coherent, Eq. ([24} (curve A), and incoherent 
(line C) light, variance for coherent (line B) and MI (curve 
D) states under scattering of incoherent light. Normal ordered 
photon-number variance for MI state under scattering of co- 
herent light is zero for all angles. Oq = 0, N = M = K = 30. 
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where the first four terms has features at angles typical 
to classical diffraction, the fourth term is also responsible 
for the doubled-frequency feature, and the last two terms 
contribute to the isotropic component. 

For the coherent state, the light scattered into a diffrac- 
tion maximum displays a very strong noise (equal to 
67V|. + Nk because {D*'^D'^) = TV^ + 6iV|- + 
Nk and {D*D) = + Nk), which is much 
stronger than the isotropic component {N]^ in highest 
order of Nk) and the features at 6i = ±7r/2 {2N^ in 
highest order of Nk) [Fig.EJb)]. In SFm state, the noise 
at maxima can be suppressed, while at other angles, in 
highest order of Nk, it is equal to that of the coherent 
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state [FiglSjc)]. In MI state, the variance (A|Dp)^ is zero 
for all angles. Conclusions about state distinguishing by 
measuring light statistics are very similar to those drown 
from the intensity and amplitude measurements, which 
have been discussed in Sec. VIIA, including scattering of 
incoherent light (see Fig. [5] and the discussion of Fig. [2]). 



In experiments, the nontrivial angle dependence of the 
noise can help in the separation of the light noise reflect- 
ing atom statistics from technical imperfections. 
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VIII. CONCLUSIONS 

We studied off-resonant collective light scattering from 
ultracold atoms trapped in an optical lattice. Measur- 
ing the light field allows to characterize the quantum 
state of atoms in a nondestructive way and in particular 
distinguish between different atomic states. The scat- 
tered light differs in intensity, quadrature variances, and 
photon statistics. A measurement of the intensity angu- 
lar distribution provides information about atom num- 
ber fluctuations in a finite lattice region, local quantum 
statistics at single sites, and pair correlations. Note that 
even local statistics can be determined by global mea- 
surements without an optical access to particular sites. 
Alternatively to angle-resolved measurements, variations 
of the mode wavelengths with respect to the wavelength 
of a trapping beam can be considered. 

Light scattering as a diagnostic tool has particular ad- 
vantages for optical lattices in contrast to scattering from 
a homogeneous BEC. Here one has a natural way to 
suppress the strong classical scattering background by 
looking at the directions of diffraction minima. In these 
directions the expectation value of the field amplitude 
vanishes while the intensity (photon number) is nonzero 
and directly reflects quantum fluctuations. Furthermore, 
in an optical lattice, the signal is sensitive not only to 
the periodic density distribution, but also to the peri- 
odic density fluctuations, giving an access to even very 
small nonlocal pair correlations. These can be obtained 
by measuring light at diffraction maxima. 

As the most striking example, we considered light scat- 
tering from a ID lattice in a transversally pumped opti- 
cal cavity as in a setup involving collective cavity cooling 
[ssl |40|. [isj . The number of photons scattered into the 
cavity is zero for the Mott insulator phase but propor- 
tional to the atom number in the superfluid phase. Both 
states have almost the same average density but different 
quantum uncertainties. So the superfluid state is a quan- 
tum superposition of different Fock states corresponding 
to all possible distributions of N atoms at M sites. Un- 
der illumination by a coherent light, various Fock states 
become entangled to scattered light states with different 
amplitudes and phases. In contrast to classical scatter- 
ing, where the atoms are described by c-number center- 
of-mass positions [1^ , for a quantum description of the 
atomic motion the light field amplitudes corresponding 
to different atomic distributions do not interfere. This 
is due to the orthogonality of Fock states forming the 
superfluid providing a sort of which way information. 

In the example configuration, the cavity-field ampli- 
tude is determined by the atom number operators fii = 
blbi. Hence, the expectation value of the field amplitude 



is sensitive to the average density only. In contrast, the 
intracavity photon number reflects the second moments 
of atom number operators (e.g. density-density correla- 
tions), while photon statistics reflects the forth moments. 

Let's emphasize that other physical systems are pos- 
sible, where the light amplitude depends on the matter- 
field amplitudes bi, while the intensity is sensitive to den- 
sity operators fii . The latter situation is typical to config- 
urations, where two or more atomic subsystems exist and 
can interact with each other, in particular, through light 
fields. The examples are matter-wave superradiance and 
amplification, where two or more momentum states of 
cold atoms were observed [ll|, [l^, 113, Ell , and interaction 
between two BECs with different internal [H, [13, HH [111 
or motional [5^ atomic states. In the framework of our 
paper, matter-field amplitudes bi can also contribute to 
light amplitudes a; , if the tunneling between lattice sites 
is important, which we have considered in the general 
model in Sec. II. In the rest of the paper, the lattice was 
assumed deep leading to negligible tunneling. 

In general, a variety of optical effects can be sensitive 
to a quantum state of an ultracold matter, if the density 
operators enter measurable quantities nonlinearly, such 
that the expectation value of those quantities cannot be 
simply expressed through expectation values of density 
operators. For instance, the x*-'^^ nonlinearity [54| and 
refractive index of a gas, where nonlocal field effects are 
important [s^ . were shown to depend on atom statistics. 
In this paper, we focused on such nonlinear quantities as 
intensity, quadrature variances, and photon statistics of 
scattered light [5^. Phase-sensitive and spectral charac- 
teristics mentioned in Sec. IV. D reflect the dependence 
of the dispersion of a medium on the quantum state of 
matter [3J,[53|- Moreover, such dispersion effects (e.g. 
cavity-mode shift) will reflect atom statistics not only in 
light intensity, but even in light amplitudes (a/). 

So far we have neglected the dynamic back action of 
scattered field on atoms. This can be well justified in a 
deep lattice where the momentum transfer is by far not 
enough to change the atomic vibrational state as long as 
not too many photons are scattered. Even without en- 
ergy transfer the information one gets from the light will 
induce measurement back-action. This should have in- 
triguing consequences for multiple consecutive measure- 
ments on the light scattered from optical lattices. 
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